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c^ ■ Abstract 

^, 

Cr. We investigate single-particle excitations and strong-coupling effects in the BCS-BEC crossover 

^ ' regime of a superfluid Fermi gas. Including phase and amplitude fluctuations of the superfluid order 

parameter within a T-matrix theory, we calculate the superfluid density of states (DOS), as well 
as single-particle spectral weight, over the entire BCS-BEC crossover region below the superfluid 
transition temperature T^- We clarify how the pseudogap in the normal state evolves into the 
,.^ , superfluid gap, as one passes through T^- While the pseudogap in DOS continuously evolves into 

Z^ ' the superfluid gap in the weak-coupling BCS regime, the superfluid gap in the crossover region 

00 

^f^ . is shown to appear in DOS after the pseudogap disappears below T^. In the phase diagram with 

^^ \ respect to the temperature and interaction strength, we determine the region where strong pairing 



fluctuations dominate over single-particle properties of the system. Our results would be useful for 
the study of strong-coupling phenomena in the BCS-BEC crossover regime of a superfluid Fermi 

gas. 
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I. INTRODUCTION 



Strong correlation between particles is one of the most important key issues in condensed 
matter physics. The recently developed ultracold Fermi gases offer unique opportunities 
to study this important topic in a controlled manner, by maximally using highly tunable 
physical parameters. Indeed, the crossover from the Bardeen-Cooper-Schrieffer (BCS) type 
superfiuid to the Bose-Einstein condensation (BEC) of tightly bound molecules [l|-y] has been 
realized in ^"Kp] and ^Li Fermi gases(6|-[8|, using a tunable pairing interaction associated 
with a Feshbach resonance 9l-ll4| . This BCS-BEC crossover demonstrates the usefulness of 
an ultracold Fermi gas as a quantum simulator for strongly correlated fermion systems. In 
particular, this system is expected to be useful for the study of high-Tc cuprates[l5|. 

The recent momentum-resolved photoemission-type spectroscopy developed by JILA 
technique to probe microscopic properties of a cold Fermi gas in the 



group 



16|] is a powerfu 



crossover region 



17 



18| . This experiment is an analo gue of the angle- resolved photoemission 
spectroscopy (ARPES) in condensed matter physics [l9l|. In the photoemission-type experi- 
ment developed by JILA group [16], atoms are transferred to the third empty atomic state 
by rf-pulse. Using this, one can directly measure the single-particle spectral weight (SW), 
as well as the occupied density of states (DOS). As a remarkable experimental result, the 
pseudogap has been observed in ^°K Fermi gas|20|. The back-bending curve of single-particle 
dispersion has been observed as a characteristic signature of pseudogap phenomenon 171. |20|. 



The pseudogap has 



of high-Tc cuprates 



15 



jeen considered as a crucial key issue in the underdoped regime 



19 



21| . So far, various mechanisms have been proposed to ex 



plain this phenomenon, such as preformed Cooper pairs 22|-|30| , antiferromagnetic spin 
fluctuations 281] , localization of Cooper pairs 2l|, |3l|, and hidden order [32]. However, the 



complete understanding has not been obtained yet, because of the complexity of high-Tc 
cuprates. In contrast, the origin of the pseudogap observed in a ^"^K Fermi gas 16|, l20|] is 
well understood. Namely, preformed pairs associated with strong pairing fluctuations are 
responsible for this phenomenon 33|] . Thus, cold Fermi gases are very suitable for the study 
of preformed pair scenario proposed in high-Tc cuprates. 

Recently, the pseudogap phenomenon above Tc has been theoretically addressed in the 



literature of cold Fermi gas 17|, 



18 



33 



-|37|. It has been shown that a 



emerges in the single-particle excitation spectra in the pseudogap regime 



gap 



33 



ike structure 



35J . It has been 



also pointed out the existence of two pseudogap temperatures T* and T** [33| : While a dip 
structure appears in DOS below T*, a double-peak structure and back-bending dispersion 
are seen in SW below T**. 

In this paper, we investigate single-particle excitations and effects of strong pairing fluctu- 
ations in the BCS-BEC crossover regime of a superfiuid Fermi gas. Extending our previous 



paper for the pseudogap phenomenon above Tc[33| to the superfiuid phase below Tc, we 
calculate DOS within a T-matrix theory. We clarify how the pseudogap above T^ evolves 
into the superfiuid gap below T^. While the evolution is continuous in the weak-coupling 
BCS regime, the superfiuid gap is shown to appear after the pseudogap disappears below 
Tc in the crossover region. We also identify the region where pairing fluctuations dominate 
over single-particle properties in the phase diagram with respect to the temperature and in- 
teraction strength. Recently, strong-coupling effects on SW has been discussed below TcJ38|. 
In this paper, we also treat this quantity to examine how the pseudogap in DOS is related 
to SW affected by pairing fluctuations. 

This paper is organized as follows. In Sec. [TTl we explain our formulation based on a 
T-matrix theory. In Sec. IIIIl we present our numerical results for the superfiuid DOS, as 
well as SW, to discuss strong-coupling effects on these quantities. In Sec. [TVt we present the 
phase diagram of a superfiuid Fermi gas to clarify the region where pairing fluctuations are 
crucial for single-particle excitations. Throughout this paper, we set ^ = A;^ = 1, and the 
system volume V is taken to be unity. 

II. FORMALISM 

We consider a three-dimensional Fermi gas, consisting of two atomic hyperfine states 
described by pseudospin a =t, i- We assume that the two hyperfine states are equally 
populated, and ignore effects of a harmonic trap, for simplicity. We also assume a broad 
Feshbach resonance as the origin of tunable pairing interaction. In this case, it is known 
that the detailed Feshbach mechanism is not crucial for the study of interesting BCS-BEC 
crossover physics, so that we can safely use the ordinary single-channel BCS model, given 
by 

H = 2_^^pCpaCpa -U 2_^ S+'7/2t'^-P+q/24.'^-P'+'7/24V+q/2t- (l) 

p,a p,p',q 



Here, cj,^ is the creation operator of a Fermi atom with pseudospin a =t, i and the kinetic 

2 

energy ,^p = ep — /i = ^— /i measured from the chemical potential /i (where m is an atomic 
mass). — f/(< 0) is a tunable pairing interaction associated with a Feshbach resonance. In 
cold atom physics, this pairing interaction is conveniently measured in terms of the s-wave 
scattering length a^, which is related to U as [39'], 

47ra, ^ U . 

where Uc is a high-energy cutoff. In this scale, the weak-coupling BCS regime and the 
strong-coupling BEC regime are, respectively, given by (fepfls)^^ ^ — 1 and (kpas)^^ > 1, 
where kp is the Fermi momentum. The region —1 < (/cpc^s) ^ ^ 1 is called the crossover 
region. 

To consider fluctuations in the Cooper channel below T^, it is convenient to rewrite Eq. 
([T]) into the form consisting of the mean-field part and fluctuation contribution. Introducing 
the Nambu field, 

%=\T 1, (3) 



we have|40l| 



H = y,Hl^pr3 - Ari]M/p - [/ V p+(q)p-(-q)- (4) 



V Q 



Here, Tj (j = 1,2,3) are Pauh matrices acting on particle-hole space. p±{q) = [piiq) ± 
ip2{q)]/2 are the generalized density operators, where Pj{q) = Xlp^ + /2''j'^p-q/2- In E<1- 
(jl]), the first term is the mean- field Hamiltonian, where the superfluid order parameter 
A = U"^ {c^piCp-^) is taken to be real and proportional to the ri-component. In this 



choice, pi{q) and P2{q) physical 



i. 



describe the amplitude and phase fluctuations of the 



411]. Namely, the last term in Eq. (jl]) describes effects of 



order parameter, respectively [40, 
pairing fluctuations. 

In this paper, we take into account the last term in Eq. (jlj) within the T-matrix 
approximation 38] . For this purpose, we introduce the 2 x 2-matrix single-particle ther- 



mal Green's function, given by 

^^^'"^"^ = GO(.u;„,)-i-S,(.c.„)- ^^^ 

Here, Gpi^iUn)^^ = iwn— ^pTs+Ari is the mean-field Green's function, where Un is the fermion 
Matsubara frequency. The 2 x 2-matrix self-energy T,p{iun) describes fluctuation corrections. 



(a) J] = 



(b) 'N/X/V 




FIG. 1: Fluctuation contributions to the self-energy Sp(iw„) in the T-matrix approximation em- 
ployed in this paper, (a) Self-energy correction, (b) Particle-particle scattering matrix Tq{ii/n)- 
The solid and dashed lines represent the mean- field Green's function Gp{iuJn) and pairing inter- 
action —U, respectively. The bubble diagrams represent the zeroth order correlation functions 
^(f i'i'^n) (where s = it), describing pairing fluctuations below Tc- The solid circles are Pauli 
matrices r=. 



Within the T-matrix theory, it is diagraminatically given by Fig. [TJ38]. Summing up the 
diagrams in Fig. [H we obtain 

q,iyn s,s'=± 

where /3 = 1/T is the inverse temperature, z/^ is the boson Matsubara frequency, and 
T± = [ti ±iT2]/2. The particle-particle scattering matrix F^'^ {ii^n) is given by 



-u 



Tq {iJ^n) rq + ('iz/„ 



1 + f/ 



n+-(2z/„) n++(zz/„; 



Fig (2z/„) U^ + {ii^r, 



(7) 



Here, 11^'' (iun) is the lowest-order of the following correlation function: 



nf(2i.„)= / dT{TApsiq,T)pAq,0)})e'''-\ 



Evaluating Eq. ([8]) within the zeroth order with respect to the last term in Eq. 



n^^ {tiy„ 



- ^ Tr \TsGl^g/^{iujn + iUn)rs>Gl_g/2{iuJn 



(3 



(8) 
we have 

(9) 



P,i^n 



Executing the a;„- summation in Eq. (|9]), we obtain 



nj+(i!'„) 



lEE 



sA2 



E.. 



p+q/2 



+ sE, 



p-q/2 



4 ^^ ^ Ep+g/2Ep-q/2 Vl + {Ep+q/2 + sEp_q/2Y 

tanh ( -Ep+q/2 j + stanh ( -Ep_q/2 



(10) 



n+-(2z/„ 



zEE 



s=±l p 



1 + S 



ip+q/2^p-q/2 
Ep+q/2Ep-q/2j *^n ~ {Ep+q/2 + sEp_q/2) 



+ 1 



^p+q/'i 



E, 



p+q/2 



1 



zi/„ 



Ep-q/2j T^l + (-£^p+<7/2 + sEp_q/2] 



tanh ( -Ep+q/2 j + s tanh ( -^p-q/2 



(11) 



where -Ep = tV^p + A^ is the Bogohubov single-particle excitation spectrum. The other 
components are given by n„^(iz/„) = n++(2i/„), and n~+(2z/„) = n+~(— «//„). 



DOS is obtained from the analytic continued Green's function, as 

p 

where SW A[p,u) has the form 

A{p,u) = lm.Gp{iijJn -^ u+ = u + i6)\u. 

n 

The analytic continued self-energy Sp(u; -^ uj+) involved in Gpiiu — t- uj+) is given by 
T,p{iu}n -^ W+)|ii = Shf 

kp+q 



(12) 



(13) 



'LpiitOn 



W+j|l2 



where z+ = z + i6, and 



q,s=±l 

-E 

qr,s=±l 

Shf = 



^ ^^UBizl+JisEp^ 

Z- sE^+a + W 



1-s- 



E. 



p+q 



Im [r+-(^+)] ,(14) 



^^ nBJz) + f{sEp+q) 



A 



Z — sEp+q + UJ+ Ep+q 



Im [r++(^+)] 



u 

2" 



-^E 



E„ 



(15) 



(16) 



is the Hartree self-energy. The other components are given by Sp(a;+)|22 = — Ep(— a;+)|ii, 
and Sp(a;+)|2i = T,p{u+)\i2. In Eqs. ( 1T^ and (TT5|) . we have carried out the Matsubara 
frequency summation by using the spectral representation of T^J^ {i^q)i given by 

q {i^n -^ Z+)\ 

(17) 



q 

r++ 



Wr, 



1 /-"^ Im [r^ 

-U-- I dz- ^ "= 



TT 



zz/„. - 2; 



«z/„ 



1 /-^ Im ^^^ 



TT 



Uz^« -> 



IVr,. - Z 




FIG. 2: (Color online) Calculated superfluid order parameter A (a) and Fermi chemical potential 
/i (b) in the BCS-BEC crossover normalized by the Fermi temperature Tp. We use these results 
in calculating DOS in Sec. III. The upper and lower insets show Tc and /i(Tc), respectively. (Tc is 
also shown in panel (b) as the dashed line.) The first-order behavior seen in the crossover region 
{{kFas)~^ ^ —0.5) is an artifact of the T-matrix approximation we are using in this paper. 



We actually calculate Eq. flT^ after determining the superfluid order parameter A and 
Fermi chemical potential /j, below Tc. In the present T-matrix theory, they are obtained by 

solving the gap equation, 

. 1 p. 

(18) 



l = f/y^tanh"P 



Ep 
2T' 



together with the equation for the number of Fermi atoms, 



^=«E^p( 



/3 



tUJn)\lie 



iSuin 



(19) 



P,(^n 



This framework is a natural extension of the Gaussian fluctuation theory developed by 



Nozieres and Schmitt-Rink (NSR)[3|], where the self-energy correction Sp(iu;„) is taken into 
account up to the first order. This T-matrix theory can properly describe the BCS-BEC 
crossover behaviors of T^ and /iJ22l. |33| . (See the insets in Fig. [2J) In addition, this theory 
is consistent with the Goldstone's theorem, in the sense that the particle-particle scattering 
matrix V^^iiVn) in Eq- O has a pole at g = z/„ = 0. Indeed, the condition that Eq. ([7]) 
has a pole at q = z/„ = gives 

[1 + f/[n++ (0) + n+=-o(o)]] [1 - f/[n++ (0) - n+jo(o)]] = o. (20) 

One finds from Eqs. (TTOj) and (ITT!) that the factor 1— ?7[nlJQ(0)— nlZg(O)] vanishes identi cally 



when the gap equation flTSj) is satisfied. We also note that, as in the NSR theory 42], 
the present T-matrix theory also shows the first-order phase transition in the BCS-BEC 
crossover region (See Fig. [21), which is, however, an artifact of the theory. To overcome this 
problem, one needs to include many-body scattering effect between molecules in a consistent 



manner 



431]. Although this is an important problem in the BCS-BEC crossover physics, in 
this paper, we leave it as a future problem and simply use A and ^ in Fig. |2] to examine 
strong-coupling effects on single-particle excitations below T^. 

III. SINGLE-PARTICLE PROPERTIES IN THE SUPERFLUID PHASE 

Figure [3] shows DOS at T^ (upper panels), as well as at T = (lower panels), in the BCS- 
BEC crossover. In the upper panels, the pseudogap structure can be seen around uj = 0|33|. 
This structure becomes more remarkable with increasing the interaction strength, reflecting 
the enhancement of pairing fluctuations. On the other hand, since thermal fluctuations are 
absent at T = 0, the well-known superfluid excitation gap associated with the superfluid 
order parameter A appears in the lower panels in Fig. [31 The goal of this section is to show 
how the pseudogap at T^ evolves into the superfluid gap below T^. 

In considering this problem, we first note the following two key issues which can be 
seen in Fig. [3l The first one is that the size difference between the pseudogap EpQ at 
Tc and the superfluid gap Eq,y at T = evaluated from Fig. [3l strongly depends on the 
interaction strength. As shown in Fig. [D^a), while E^q is smaller than E^-p in the weak- 
coupling BCS regime, the former becomes larger than the latter in the BCS-BEC crossover 
region. This implies that, while the pseudogap in the BCS regime may smoothly change 

8 
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FIG. 3: (Color online) Upper panels: Pseudogap in DOS p{uj) at Tc. Lower panels: DOS at T = 0. 
In this figure, the left and right panels show the results in the BCS side {{kpas)^^ < 0) and in 
the BEC side {{kpas)^^ > 0), respectively. At T = 0, small but finite intensity at w ~ seen in 
panel (bl) is due to a small imaginary part {6 = O.OlTp) introduced to the energy in numerical 
calculations. 



into the superfluid gap below T^ the large pseudogap at Tc in the crossover region needs to 
shrink below Tc, in addition to the opening of the superfluid gap. 

The second key issue in Fig. [3] is that, although the pseudogap structure in DOS looks 
similar to the superfluid gap, the former does not have the coherence peaks at the gap edges. 
However, even in the superfluid phase below T^., the coherence peaks are known to disappear 
by strong-coupling effects. In the present case, these strong-coupling effects involve pairing 
fluctuations excited thermally and the formation of tightly bound molecules. While the 
former effects are expected only at finite temperatures, the latter may exist down to T = 
in the BEC regime. The latter effect can be easily confirmed by using the superfluid DOS 




FIG. 4: (Color online) (a) Comparison of the pseudogap size EpQ at Tc and superfluid gap size 
£'sF at T = evaluated from Fig. [3l For comparison, we also show the energy gap Eq in the 
BCS-Leggett crossover theory _4] (which equals A when ^ > 0, and equals -^/Jj^ + A^ when /x < 0). 
The panel (b) shows the behaviors of these quantities in the strong-coupling BEC regime, where 
one finds -Epc — E^f — Eq when (kpag)^^ > 3. The inset shows how to determine EpQ and E^f 
from Fig. [3l Since the pseudogap actually does not have a clear energy gap, we conveniently define 
the gap size as the half of the dip size at [p{too) + p('^i)]/2, as shown in the inset, where ujq and wi 
are the bottom energy and lower peak position, respectively. 



10 



within the mean-field theory, given by 



p{uj) 






[e{iu - A) - e{-uj - A)] 



Vuj^ - A2 + 



/^ 



u 



v^ 



u^ 



/^ 



m 



3/2 



Vw2_A2 



A2 

(/i > 0), 



2v^7r2 

X 



Vaj2_A2 



(/i < 0). 



(21) 



Vw2 _ A2 



l/^l 



In the BCS regime where /x > 0, the singularity in u/^/u^—A^ in the upper equation gives 
the diverging coherence peaks at cj = ±A. In contrast, the singularity at l^l = A in the 
lower equation is less important in the BEC regime when yU < 0, because DOS is finite only 



when 1^1 > a/A^ + /i^ > A. Since the negative /i in the BEC regime is a strong-coupling 
effect associated with the formation of tightly bound molecules|2|, yj, the suppression of the 
coherence peaks in the BEC regime at T = may be also regarded as a strong-coupling 
effect. From the above discussion, we find that the coherence peaks in DOS may be used to 
determine the region where strong-coupling effects are less important and one can discuss 
superfluid properties to some extent within the weak-coupling mean-field BCS theory. 

We now consider the superfluid DOS p{ijj) below Tc. Figures [5E show the temperature 
dependence of calculated p{u) in the BCS-BEC crossover. In the weak-coupling BCS regime 
(Fig. |5]), the pseudogap at Tc is found to smoothly change into the superfluid gap below Tc. 
One can see the growth of the coherence peaks in panels (a2) and (a3), and the BCS- 
type superfluid DOS having a clear gap structure with sharp coherence peaks is eventually 
realized far below T^ (panel (a4)). We note that a similar continuous evolution from the 
pseudogap to the superconducting gap has been observed in the underdoped regime of high- 
Tc cuprates|2ll. |46|. 

Figure [5] also shows the intensity of SW in the right panels. As discussed in our previous 
the pseudogap phenomenon in SW is not remarkable in the BCS regime. Indeed, 



paper 



a peak line corresponding to the free particle dispersion u = p^ /2m — ^ is only seen in panel 
(bl), although the pseudogap can be clearly seen in panel (al). The expected superfluid 
gap simply opens at w = below Tc, as shown in panels (b2)-(b4). Comparing these results 
with the mean-fleld expression, 

1 



A{p,u) 



1 + 



/i 



E. 



p J 



6{uj - Ep) + 



fi 



E., 



p J 



6{lo + Ep) 



(22) 
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FIG. 5: (Color online) Temperature dependence of DOS p(a;) and intensity of spectral weight 
Aijp^ijS) in the weak-coupling BCS regime {(kpas)~^ = —1). 

we find that the overall behavior of SW in the BCS regime is essentially the same as that 
in the mean-field theory. 

We obtain quite different results in the unitarity limit shown in Fig. |6l In this case, 
panel (a2) clearly shows that the superfiuid gap structure is still absent around w = even 
at T = 0.9Tc. This is because, although the superfiuid order parameter A itself is finite, 
the superfiuid gap structure in DOS is smeared out by strong pairing fluctuations at this 
temperature. However, panel (a2) also shows that pseudogap structure becomes obscure, 
indicating the suppression of pairing fluctuations (although they are still strong enough to 
smear out the superfiuid gap in DOS). 

This suppression of pseudogap below T^ can be also seen in SW. In Fig. Et^bl), we can 
see the typical pseudogap structure of SW, namely, the double-peak structure consisting 
of a positive energy (particle) and negative energy (hole) branches. At T = 0.9Tc (panel 
(b2)), this double-peak structure becomes obscure due to the appearance of finite spectral 
intensity around w = 0. However, the superfiuid gap still does not open at cj = at this 
temperature, being consistent with panel (a2). 
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FIG. 6: (Color online) Same plot as in Fig. [5] for (kpUs 



(unitarity limit). 



At lower temperatures, when the pseudogap in DOS almost disappears, Fig. |6t^a3) shows 
that a dip structure appears around u = 0. Correspondingly, SW also has a gap structure 
at cj = 0, as shown in panel (b3). These superfiuid gap structures develop at lower temper- 
atures, and they eventually reduce to the BCS-type DOS and SW far below T^, as shown in 
panels (a4) and (b4), respectively |45|. 

The above results indicate that, in the crossover region, the superfiuid gap appears in DOS 
after the pseudogap almost disappears. This is quite different from the continuous evolution 
from the pseudogap to the superfiuid gas in the weak-coupling BCS regime. Since strong 
pairing fluctuations, which is essential for the pseudogap phenomenon, must be suppressed 
to obtain the superfiuid gap in DOS and SW, the evolution from the pseudogap to superfiuid 
gap is a competing phenomenon in the crossover region. 

In the strong-coupling BEG regime, since tightly bound molecules have been already 
formed far above Tc, DOS has a clear gap structure even at Tc, as shown in Fig. El^al). 
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FIG. 7: (Color online) Same plot as in Fig. [5] for {kpas 



0.8 (BEC regime). 



Because of this clear gap structure, although the shrinkage of this gap can be seen in panel 
(a2) (which is considered to correspond to the suppression of the pseudogap discussed in 
Fig. E]), one cannot precisely determine the temperature where the superfluid gap structure 
at cj = starts to appear in Fig. [71 In this regard, we note that the chemical potential /z is 
negative in the case of Fig. [3 Thus, as will be discussed in Sec. IV, one should regard the 
system in this regime as a molecular Bose gas, rather than a Fermi gas. 

Figures El^bS) and [7l^b4) show the appearance of a sharp negative energy branch in SW 
far below T^. This means that the overall spectral structure becomes close to the mean- 
field superfluid result given by Eq. ( I22l) far below Tc (although the chemical potential 
/i remarkably deviates from the Fermi energy e-p.) For the appearance of this negative 
energy (hole) branch, we briefly note that it is absent in the BEC limit where the molecular 
formation occurs within the simple two-body physics. On the other hand, SW has both a 
positive energy (particle) and negative (hole) branches in the BCS regime, reflecting that 
Cooper pairs are many-body bound states assisted by Fermi surface. Thus, the sharp hole 
branch in panels (b3) and (b4) indicates that the many-body effect still contributes to pair 
formation to some extent even in the BEC regime at {k-pas)"^ = 0.8. 
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In the strong-coupling BEC limit, the Green's function in Eq. ([5]) reduces to 38|, |44| 



«-<-'i-;4TiTA|^- (^^' 



(We summarize the derivation of Eq. ( !23|) in the appendix.) Here, 






is the superfluid order parameter in the BEC limit, and ApQ = ^JSTxh-Q/ni^as is the pseu- 
dogap parameter |38|, where rig and h-Q represent the molecular condensate density and 
molecular non- condensate density, respectively. Equation ( 123|) shows that the single-particle 
excitation gap Eq is given by 



E.^,/,^ + (A= + A|,)^,/(^) +^. (25) 

where n = 2[ng + ne] is t^ total fermion density, and we have used the expression /i 
-l/2ma1 in the BEC limit 



Equation fl25|) means that the excitation gap in DOS becomes 
T-independent deep inside the BEC regime. Indeed, Fig.lD^b) shows that Epc — -Esf(— Eq) 
when (kpas)'^ > 3. In the extreme BEC limit (a~^ — )■ +cxd), Eq. ( I25l) reduces to the half of 
the binding energy of a two-body bound state -Ebind = ^/rnaj., as expected. 

Before ending this section, we note that, although the overall structure of DOS is very 
close to the BCS-type DOS at T = (See the lower panels in Fig.O), the gap size Esf in DOS 
at T = is smaller than the magnitude of superfluid order parameter evaluated in the BCS- 
Leggett crossover theory [2|] (which consists of the mean-field gap equation and mean-field 
number equation), as shown in Fig. HJ^a). This is because the self-energy correction T,p{iun) 
in Eq. ([5]) still affects single-particle excitations even at T = 0, although pairing fluctuations 
are suppressed far below T^. Indeed, in the present T-matrix theory, the superfluid gap in 
DOS is affected by the off-diagonal self-energy Sp(iw -^ u + i6)\i2 at cj ~ A even far below 
Tc. In addition, the present strong-coupling theory involves effects of an effective molecular 
■Meraction withuUhe Born approx,matio„. uarnely. the effective mofeeular seatteriug fength 
equals ob = 2as|39!] in the BEC regime. This effective boson-boson interaction leads to the 
quantum depletion n^ (which describes the number of non-condensate fermions at T = 
0), which decreases the condensate fraction n^ = [n — nj]/2, as well as the magnitude 
of superfluid order parameter in the BEC regime given by Eq. (|24|) . These effects are 

15 



completely ignored in the BCS-Leggett crossover theory [2j, so that the superfluid gap size 
EsF in DOS becomes smaller than A evaluated in the mean-field-type crossover theory. 
However, since these strong-coupling effects are eventually suppressed deep inside the BEG 
regime, the difference between the two theories becomes small in the EEC limit, as shown 
inFiglUb). 

IV. PHASE DIAGRAM IN THE BCS-BEC CROSSOVER 

Figure E] shows the phase diagram of a cold Fermi gas in the BCS-BEC crossover. In 
panel (a), we introduce three characteristic temperatures, T*, T*, and 2|/i(Tc)|, in order 
to conveniently identity the region where pairing fluctuations dominate over single-particle 
properties. T* is the temperature where the superfluid gap appears in DOS below Tcj47|. 
The region above T* is considered to be dominated by strong pairing fluctuations even in 
the superfluid state. T* is the so-called pseudogap temperature discussed in our previous 



paper J33|, where the pseudogap starts to emerge in DOS above Tc. In addition, we also take 
into account the fact that physical properties in the strong-coupling BEC regime are close to 
those of a molecular Bose gas, rather than a Fermi gas. Noting that the molecular binding 
energy -Ebind in this regime is deeply related to the Fermi chemical potential as -Ebind — 2|/i| 
when yU < 0, one may expect that the thermal dissociation of molecules is suppressed in 
the BEC regime when T < 2|yu| (/i < 0). Thus, it is convenient to regard the right side of 
the 2|yu(Tc)|-line in Fig. [8](a) as the molecular Bose gas regime |33f|. We briefly note that T*, 
T*, and 2|/i(Tc)|, are all crossover temperatures without being accompanied by any phase 
transition. 

Using these three characteristic temperatures in Fig. [8]^a), we flnd that the region "PC" 
in Fig. MJo), which is surrounded by T*, T*, and 2|/x(Tc)|, is the one where the pseudogap 
structure in DOS is remarkable. Thus, we conveniently call this region the pseudogap region 
(although, strictly speaking, the region below Tc is the superfluid state). In this pseudogap 
regime, strong pairing fluctuations induce a gap-like structure in DOS in both the normal 
and superfluid phases. 

Below T* ("SF" in Fig. [8](b)), instead of the disappearance of the pseudogap, the su- 
perfluid gap starts to develop in DOS, so that single-particle properties are dominated by 
superfluid gap. As one further decreases the temperature below T*, one eventually obtains 
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FIG. 8: (Color online) Phase diagram of a cold Fermi gas in the BCS-BEC crossover, (a) Char- 
acteristic temperatures introduced in this paper. T is the temperature below which the superfluid 
gap appears in DOS 



47| . 2|/u(rc)| in the BEC regime (/i < 0) gives a characteristic temperature 
below which thermal dissociation of bound molecules is suppressed. Thus, the right side of this 
line may be regarded as the region of tightly bound molecular Bose gas, rather than a Fermi gas. 



331], where the pseudogap structure starts to 



T* is the pseudogap temperature obtained in Ref. 
appear in DOS above Tc. In addition to these characteristic temperatures, we also introduce Tcp 
as the temperature at which the BCS-type coherence peaks appear in DOS. Below Tcp, the system 
becomes close to the simple weak-coupling BCS state, at least with respect to single-particle ex- 
citations, (b) Phase diagram of a cold Fermi gas. PC: pseudogap phase. NF: normal Fermi gas. 
NB: normal state molecular Bose gas[33]. SF: superfluid Fermi gas with a superfluid gap in DOS. 
MBEC: BEC of molecular bosons. We emphasize that only Tc is the phase transition temperature, 
and the others are all crossover temperatures without being accompanied by any phase transition. 
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the weak-coupling BCS-type DOS characterized by a clear excitation gap and coherence 
peaks, as discussed in Sec. III. To conveniently include this, we also introduce the charac- 
teristic temperature Tcp at which the coherence peaks appear in DOS|48| in Fig. [Hl^a). Below 
Tcp, single-particle properties are close to those in the weak-coupling mean-field BCS state. 
In the molecular Bose gas regime ("NB" and "MBEC" in Fig. [H](b)), a large single- 
particle excitation gap already exists above Tc, refiecting a large molecular binding energy 



{Eund — 2|/i| ^ l/ma^[39|]). This large binding energy suppresses single-particle excitations 
accompanied by pair breaking in the superfiuid phase below Tc, so that excitations are 
dominated by collective Bogoliubov modes, as in the case of Bose superfiuid. In this sense, 
we call the superfiuid region in the molecular Bose gas regime the BEG of molecular bosons 
("MBEC" inFigEKb)). 

V. SUMMARY 

To summarize, we have investigated single-particle excitations and strong-coupling effects 
in the BCS-BEC crossover regime of a superfiuid Fermi gas. Extending our previous work 
above Tc to the superfiuid phase below Tc, we have numerically calculated the superfiuid 
DOS, as well as SW, within the T-matrix theory. We have systematically examined how 
the pseudogap at Tc evolves into the superfiuid gap, as one decreases the temperature below 
Tc. While the evolution is continuous in the weak-coupling BCS regime, the superfiuid gap 
was shown to appear in DOS after the pseudogap is suppressed below Tc in the crossover 
regime. Using these results, we have identified the pseudogap region where strong pairing 
fiuctuations dominate over single-particle properties in the phase diagram of a cold Fermi 
gas. 

Since the observation of single-particle excitations has recently become possible in cold 
Fermi gases by photoemission-type experiment, measurements of single-particle excitation 
spectrum affected by strong pairing fluctuations discussed in this paper would be an in- 
teresting problem to understand the strong-coupling superfluid properties in the BCS-BEC 
crossover. 

In this paper, we have assumed a uniform Fermi gas, for simplicity. In a trapped system, 
it is an interesting problem how the spatial inhomogeneity affects the evolution from the 
pseudogap to the superfluid gap below Tc. Since a real Fermi gas is always trapped in a 



18 



harmonic potential, this is also an important issue in comparing experimental data with 
theoretical calculations. We will discuss this problem in a future paper. 
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Appendix A: Analytic results in strong-coupling BEC limit 

In this appendix, we p resent the outline of the derivation of Eq. f l23p . For more details, 
44| . In the BEC limit, the particle-particle scattering matrix in Eq. 



we refer to Refs. 



([7]) reduces to [381. 



38, 

44| 



'(Al) 



where E^ = \/ j—{4 — h 2|yUB|) is the Bogoliubov excitation spectrum of a molecular BEC, 



<7 V 4m *> Am 

>2 



and fiB = ~^ /4|/i| is the Bose chemical potential. The number equation in the BEC 
limit reduces to A^ = 2 (rig + ub), where n^ and ne are the molecular condensate and 
non-condensate density, respectively, given by 

< = ^y^^, (A2) 

ly^^u,,+l/^rn-J^B_ 

Using Eq. flAip . we approximate the self-energy in Eq. (|6]) to 

^ r \ Svr I nBGl{iuJn)\22 n%Gl{iu}n)\i2 \ ,.^, 

Sp(za;„) ~ ^— , (A4) 
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where we have approximately set q = Un = in Gpj^_q{iijJn + ii^n)- Using Eq. (lA4p . one 
obtains the diagonal component of the Green's function as 

Gp{iu„)\u = 1 — r2-- (A5) 

SnnV, ^n / ■ m 

lUn + 4p 5— G (2W„) 11 

Expanding the denominator in Eq. (lASp up to O(A^), we obtain Eq. (l23l) . 
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